A mathematical model of DNA degradation: Possible 
role of magnetic nanoparticles 



V. N. Binhi 

General Physics Institute, Russian Academy of Sciences 
38 Vavilova St., Moscow, 119991, Russian Federation. 
E: binhi@kapella.gpi.ru 



Abstract 

A mathematical model of genome degradation is proposed that takes into account 
a variable rate of mutation and increasing number of cells in a developing human 
organism. The model explains known properties of cancer development, in particu- 
lar, a synergism between different mutagens and an increased probability of cancer 
in the early years of life. An iteration equation is suggested that uses only a few 
model parameters and describes basic regularities observed in cancer onset. In the 
model context, relatively small chronic variations in the intracellular content of free 
radicals may markedly affect the probability of a cell to become a cancer cell. On 
the other hand, magnetic nanoparticles are shown to be an endogenous source of 
chronic magnetic exposure that increases the local concentration of free radicals. 
An enhanced level of leukaemia in early childhood is assumed to originate from 
magnetic nanoparticles located in hematopoietic stem cells. 

Key words: genome degradation, leukaemia, magnetosome, static magnetic field, 
hematopoietic stem cell, mathematical model, superparamagnetic nanoparticle, can- 
cer risk factor, radical pair mechanism 

1 Introduction 

A number of studies have demonstrated that magnetic nanoparticles found in living tissues, 
including the human brain, are involved in biological reactions to magnetic fielcS exposures 
[59] . However, the possible role of magnetic nanoparticles in molecular processes underlying 
cancer development has not yet been discussed. 

It has been repeatedly demonstrated that organisms can biochemically precipitate minerals 
[7i] including magnetic minerals like magnetite Fe304 [THl [12l |32], maghemite 7-Fe203 [98] . 
and greigite Fe3S4 [7^. The presence of magnetic nanoparticles is well documented in many 
organisms [6] . They are also found in the human brain [611 Ell El] and other human tissues 



Membrane-enclosed crystals of magnetite are often called magnetosomes [1]. Efforts were 
undertaken to identify genes required for magnetosome synthesis and arrangement. Some of 
the genes have been suggested to play a specific role [HI HHl [HI [THl [331 El]- 

Magnetosomes are often assumed to underlie the observable biological effects from expo- 
sure to weak MFs |66j. The energy of a 100-nm magnetosome in the geomagnetic field is 
~ 24 /cbT. When exposed to an additional variable magnetic field /i, the energy varies in 



-'^Abbreviations: MF (magnetic field), HSC (hematopoietic stem cell), DNA (deoxyribonucleic acid), SD 
(single domain), SP (superparamagnetic), RP (radical pair). 
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about (/i//fgco)24 /cbT. If these changes exceed thermal fluctuations ~ A;bT/2, they can cause 
a biological response. This sets a natural constraint on the MF magnitude capable of affecting 
a biophysical or biochemical system appreciably: h > 1-2 /iT. Recent studies have shown that 
nonlinear stochastic dynamics of a magnetosome may be a basis for explaining biological effects 
from MFs of as low as 200 nT P], and of biological magnetic navigation with a high accuracy 

The biological role of magnetic nanoparticles in higher organisms is not completely under- 
stood. The question is whether the consequences of their presence in biological tissues are 
limited only to magnetic orientation and possibly navigation by migrants, or can magnetic 
minerals take part in destructive processes in cells? 

There are some indications that biological cells are not indifferent to magnetic nanoparti- 
cles. As was first noted in [8j , MFs produced by magnetic nanoparticles are orders of magnitude 
greater than the geomagnetic field, and this can be an important endogenous source of chronic 
magnetic exposure facilitating free radical formation around the particles. Possible medical 
implications of human contamination by magnetic compounds are discussed in [HS]- Experi- 
mentally, ultrafine particles 12-14 nm in size were shown to be internalized by human monocyte 
cells and significantly increase, by 40-45%, the release of free radicals |113j . The severity of 
neurodegenerative diseases has been found to correlate with the amount of magnetite in the 
human brain [SlHHlETj, which is explained by excessive toxic ferrous ions in brain. Formation of 
weakly charged ferric oxide nanoparticles in DNA complexes was discussed in [H^ with regard 
to observation of broad ESR spectra in DNA preparations. 

It is assumed in the current article that some higher organisms have an abnormally in- 
creased content of magnetic nanoparticles in hematopoietic stem cells (HSCs). The content is 
low enough not to be detected by ordinary magnetic methods, since it is masked by the para- 
magnetic contribution of blood and the superparamagnetic contribution of ferritin. However, 
it is high enough to contribute to DNA lesions by chronically producing an additional amount 
of free radicals. HSCs are thought to be one of the primary places where acquired mutations 
can accumulate [53] . 

The aim of this study is, first, to show within a mathematical model that even a small 
chronic change in the intracellular concentration of genotoxic agents such as free radicals can 
markedly enhance the rate of cancer incidence in the early stage of an organism's development 
and, second, to show that chronic changes in the free radical concentration can be induced by 
intracellular superparamagnetic nanoparticles. 

Many mathematical models have been proposed to describe different aspects of cancer de- 
velopment: the evolutionary multistage model of cell proliferation in a changing environment 
[33] . the multiscale model of key genes, cellular kinetics, and tissue dynamics [99j, the model 
involving the dynamics of gene inactivation with genetic instability [HHl EH] , and others recently 
reviewed in [63, [16]. Usually, the rate of mutation is considered to be a parameter. However, 
the rate of mutation may be a nontrivial function of the rate of DNA lesions, depending on cell 
cycle duration and the time needed to repair the lesions. 

In this paper, a mathematical model is proposed that takes into account a variable rate of 
mutation. The model explains many known properties of cancer development: in particular, a 
synergism of chemical and radiation mutagens, and an increased probability of cancer in early 
years. It also provides an unobvious correlative link between the risk of getting some cancers 
in middle age and the time when acute exposure to an exogenous genotoxic factor took place. 
Numerical estimates are proposed for possible increases in the rate of cancer incidence from 
acute and relatively small chronic changes in intracellular free radical content. 

Some recent experiments on the biological effects of static magnetic fields are summarized, 
average proper magnetic fields generated by superparamagnetic nanoparticles are calculated. 
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and an idealized free radical reaction is analyzed to show that a magnetic nanoparticle can 
increase the rate of free radical formation in a cell by a few percent, which is enough to cause 
marked consequences in the rate of DNA degradation. 

2 DNA degradation model 

DNA molecule is a sequence of its relatively short parts, or genes, each of which encodes a 
specialized protein. Since a DNA is a very long chain, it can be damaged in one or several 
places, and often this happens without loss of its ability to replicate. If damages, or mutations, 
are not repaired to the beginning of next cell division they may become permanent mutations 
inherited in daughter cells. Mutated genes may cause the production of abnormal proteins 
leading to diseases, in particular to cancers. 

A cancer develops in a process that involves mutations most likely in a few from about 
a hundred known cancer-related genes [116] . These key genes, called oncogenes and tumor 
suppressor genes, take part in providing a dynamical balance between growth and death of 
cells of a given type. 

Various agents are known to induce mutations in living cells: many chemicals, aggressive 
free radicals, radioactive and ultraviolet radiations. Spontaneous mutations also happen. 

Normally, gene mutations are repaired by the DNA repair system. It includes a number of 
pathways through which an absolute majority of different damages in DNA may be recognized 
and repaired. The number of unrepaired mutations per year is only about a dozen in a cell. 
DNA repair system is a very perfect system that protects the body from diseases and untimely 
senescence. However, there are special genes encoding proteins, which regulate the DNA repair 
system, and these genes may also be corrupted by mutations. Then the rate of unrepaired 
mutations dramatically increases and a cancer becomes almost inevitable. 

Most cancer researchers are focused on mutations in a small number of cancer-related genes 
as the primary damages which cause the subsequent malignant transformation of normal cells 
to cancer cells. There are other theories, in which some chromosomal abnormalities, like uneven 
number of chromosomes, are a root cause of every cancer [88l [37j . However, all agree that gene 
mutations are required. 

Evidently, mutations may accumulate only in proliferating cells. There are limited number 
of such cell types. HSCs play significant role in the human body. They are pluripotent cells: 
all other blood cells are made from adult stem cells in the bone marrow [70] . Like other blood 
cells, new HSCs are also originate from adult stem cells and gradually replace them in their 
functioning. It is the reason why HSCs are thought to be one of the primary places where 
acquired mutations can accumulate. The mutations are inherited in daughter stem cells, and 
this continuous self-renewal is accompanied by a consecutive DNA degradation. The gradual 
DNA degradation in HSCs is associated with the increased probability of leukaemias and, 
along with genome degradation in other cells, with natural ageing. If a hematopoietic stem cell 
became cancerous, it produces blood cells that further proliferate out of control and result for 
example in chronic or acute myeloleukaemias. 

Below we will consider a mathematical model that illustrates main features of DNA degra- 
dation, to finally show why magnetic nanoparticles may be thought of as a cancer risk factor. 
In this model a cell experiences a set of n consecutive renewals, which fix mutative damages 
accumulated during a given cell cycle. The probability P„ of the cell to be a cancer cell at the 
end of cycle is a basic variable of the model that features the current state of the cell. 

We will study first the process of the DNA degradation within the time interval T of a single 
cell cycle. Because only about a hundred from 35000 genes in human genome are associated 
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with cancers, and only a few of them, key genes, may relate to leukaemia, we consider one of 
these few genes and evaluate the probability p that this gene is broken to the end of a cycle, if 
it was intact before the cycle. We will assume that a gene is a set of k bits, and any damaged 
bit means the gene is damaged. 

Let mutative DNA lesions be a Poisson process defined on a statistical ensemble of virtual 
identical bits. The probability for a lesion not to occur within the time interval [0, t] is exp{—at), 
where a is the mean rate of lesions, both spontaneous and induced. 

Let t be the moment of a DNA lesion and r be the time interval needed for the cell to repair 
any damage. In this idealized scheme, lesion is a random event and reparation is a deterministic 
process, which is necessary to occur if the cell has enough time before genome duplication, or 
mitosis, starts. Consequently, if t < T — t, reparation occurs and no mutation appears in the 
gene. Oppositely, if t > T — t, the mutation will take place. 

The probability that a lesion will occur in the range [t, t + dt] is a product of two events: (i) 
the lesion had not occurred to the moment t and (ii) it has occurred within the time interval 
dt. It follows that the probability of the lesion occurring within [t, t + dt] equals exp(— Q;t)Q;dt. 
So the probability that a bit will be damaged in the range [T — r, T], i.e. the probability of 
mutation. 



Of course, the sum of the probabilities of the lesion to occur (i) before T — r, (ii) after T, and 
(iii) p is equal to unity. Here it is assumed that more than one lesion in the same bit per cell 
cycle is an unlikely event. 

Next, we define that if the reparation time r is greater than the cell cycle duration T, the 
probability p no longer depends on r and remains at that level where r = T, i.e. 1 — exp(— aT). 



This function is plotted in the Fig. [H 

The probability that a gene of k bits is undamaged to the end of a cell cycle is (l—p)^, and 
that it is damaged 1 — (1 — pY . 

Cells carrying the damaged genes cannot work properly. Most of such cells undergo apopto- 
sis or the immune system will detect and kill them. However, a very small part of such mutant 
cells retain their functions and continue dividing. Let u be the probability that a cell carrying 
a damaged gene will enter the next cell cycle. Then, an acquired mutation appear with the 
probability p = u[l — {1 — p)^]. The probability of no acquired mutation is g = 1 — p to the end 
of a cell cycle, and Qn = q'^Qn-i to the end of cycle n in a population of cells. Actually, 
both q and A^ depend on n, and so Qn = qn"Qn-i- Then the probability of the population to 
have an acquired mutation of a key gene in a cell to cycle n is 



Pn = 1 - = 1 - (1 - Pn-l)il-Pnf- = 1 - (1 - P„_i) {l - u[l - (1 . (2) 



This quantity, which is associated with the probability of cancer, is a subject of the subsequent 
study. 

As was said, the DNA repair system may be corrupted by mutations. In the model, there 
is the only parameter that describes the efficiency of the repair system, the reparation time r. 
The dependence of reparation efficiency on mutations means that r is not a parameter but a 
variable r^. We assume as a model approach that the repair system depends on mutations to 
different key genes in the same manner, irrespective of gene specificity. This enables to farther 




Then 




(1) 
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Figure 1: The probability that a bit of a gene is damaged to the 
end of a cell cycle. To the left of T/r = 1 the proba- 
bility increases that the cell has enough time to repair 
a damage; to the right of T/r = 1 lesions have no time 
to occur, a is the mean rate of lesions. 

simplify the model: the only key gene will be studied, and the probability of it to be damaged 
reflects that of the cell genome as a whole. In this case, the variable r„ depends on P„. 

By definition, the value P„ increases monotonously between and 1. At P„ = 0, the 
reparation time should be equal to its initial value r(0) = u. At Pn I, Tn should tend to 
a maximum, for which T is a suitable value, as was discussed above. In addition, r„ should 
depend on the value of P at the previous cell cycle. As a reasonable approximation, we take a 
linear dependence r„ = + P„_i(T — z/). Substituting this in ([T]) and using the designations as 
follows 

a = aT , b = ah',c = a — b>0, (3) 

we write, for the case t < T, which takes place almost always except a special situation 
considered later, 

Pn = b + cPn-l + 0{b)+0{c), (4) 

where a <^ 1, 6 <^ 1. In what follows, the mathematical indications of the accuracy of approx- 
imated equalities are omitted for convenience. Given this, we use the sign of exact equality 
keeping in mind that only terms linear in the small parameters a, b, u, and others are saved. 
Where these parameters appear as products like akN, which are not necessarily small, it will 
be specially indicated. 

Substituting the latter equation in gives an iteration equation 

P„ = 1 - (1 - P„_i) {1 - n[l - (1 - 6 - cPn-i)']f" . (5) 

In order to study the behavior of the probabilities at large times, we assume that the number 
of cells Nn has reached its stationary value N. The results of computation of equation are 
shown in Fig. [2] along with epidemiological data on cancer rates in the USA. Calculations were 
made with a = 10~^ and k = 2x 10^. For "leukaemias" case other parameters were = 5 x 10^, 



5 




Time t, days 



Figure 2: Probabilities P that a key gene mutation appears in 
a population of cells to the time t calculated from 
equation ^ with = 5 x 10'^, curve 1, and = 5 x 10®, 
curve 2. Points: cancer rates in the USA, all races, male 
and female, adapted from SEER 13 Registries for 1994- 
2003, http:/ /seer, cancer .gov. 



T = 10, z/ = 2 X 10-^ M = 1.6 X 10-^ For "all sites" case, A^ = 5 x 10^ T = 100, v = 2x 10-^ 
and u = 2.2 x 10~^, all times in days and time variable t = nT. Parameter v controls the 
slope of the curves in the inflection area. Parameter m, related to the immune and apoptotic 
processes, and z/ were fitting parameters, and the product au rather than a and u separately 
was significant for calculations. Numerical values found from this fit are discussed below and 
used in what follows as reference values around which the quantities of individual organisms 
may vary. 

As is seen, there are three distinct areas of curves' behavior. It is slow quasistationary raise, 
an inflection area, and fast approaching to unity. On the whole, such behavior is associated 
with a slow-growing probability of cancer until the inevitable genomic degeneration against the 
background of senescence-related loss in the quality of genetic information. 

The time position of curves completely depends on the value of P taken as initial value 
for the iteration process at large times, with all other parameters fixed. When Pini is 10 times 
greater, probability of cancer starts to grow significantly by 30 years earlier. Fig. [3], i.e. about 
10-year shift per duplication of Pini. Formation of this value of P at times about several years 
occurs in early stage of population development, along with acquired mutations. It is interesting 
to find an analytical solution of equation ^ in the range of relatively small values of n, since 
it is the interval where genetic mutations rapidly accumulate. 

In this case P„ <^ 1 and equation ([5]) may be simplified. Based on digital values used for the 
parameters in the above calculations, one can see that kb^l and (1-6-cP)^ = 1-kb-kcP, 
then 

P„ = 1 - (1 - P„_i) {1 - Uk[b + cPn-l]f" . 
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Figure 3: Shift of probability curve due to variation of the initial 
value of P in iteration equation (l5l). Points: cancer rates 



in the USA, adapted from http://seer.cancer.gov 



As well, ukNc <^ 1 that yields the next equation 

Pn = Pn-l+ ukbNn, 

which approximates ([5]) with an accuracy better than 1%. Since the increment of P„ at each 
iteration step is very small, it is convenient to reduce this iteration equation to a differential 
equation. Formally replacing n with a continuous variable 6 and Pn with P{6), we write 
dP = ukbNdO. Since b = au and all the parameters except k may depend on 6, the solution 
reads ^ 

Pi9) = k [ auuNde, (6) 
Jo 

where for simplicity 6 is used also as an integration variable. It is an important relation 
demonstrating the basic regularities observable in cancer development: 

(i) Genomic degeneration, or accumulation of acquired mutations, has a character of perfect 
integration. Consequently, the effects of exposures to even a short-term mutative factors will 
never abate and will necessarily contribute to the overall degradation of the organism. 

(ii) There are different channels by which the rate of genetic degradation might be influenced. 
It is a (rate of lesions), z/ (rate of DNA repair), u (immune responses and apoptosis), and 
(number of cells and cell cycle duration). 

(iii) External factors affecting genome through the channels of influence and doing it even 
relatively weakly but chronically may significantly expedite the rapid growth of cancer risk. 

(iv) Contributions of small changes in different channels are summed up arithmetically. 
However, there is a synergistic effect of relatively strong simultaneous influences: the contribu- 
tions are multiplied in this case. For example, short-term 10-fold simultaneous growth of any 
two parameters gives rise to 100 times greater degradation rate. Probably, free radicals are so 
dangerous as they affect simultaneously many channels. 

The time variable 6 in equation ([H]) is a dimensionless variable, a time measured in the 
units of cell cycle duration T. Generally, T depends on time. It is well known that many cell 
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types divide at a varying rate so that T being minimal in the beginning then increases until it 
reaches a stationary value or even tends to infinity when cells stop dividing. This introduces 
an interesting feature concerned with the possibility that T goes through the value of r where 
the probability p that a gene is damaged to the end of a cell cycle is maximal, Fig. [H 

In this case, the approximation r < T, used in (jl]), does not work and Pn, which now does 
not depend on P, in the initial phase where P <^ 1, should be written as 



p 



aT, u > T 



which follows from ([T]). Substituting it into ([2]) with ukNb ^ 1 we obtain instead of (E]) the 
equation 



P(e) = k [ p{a,u,T)uNde, 
Jo 



where all the variables are functions of 6. We will change variable 6, i.e. dimensionless time 
expressed in T units, to time t. As the cell cycle time T changes with time t, the number of 
the cell cycles n and as well 9 completed to a moment t will not be directly proportional to t: 



nt) 

Consequently, since the function under the following integral does not depend on t explicitly, 

P = k puNde = k / puN—dt, 
Jo Jo 

or 

Pw^tjr' ^'°-y)"^ d.. (7) 

where all the variables are now functions of physical time t. To analyze this equation, some 
ideas on the character of functions N{t), T(t), and u(t) should be formulated. 

There are no age-specific data on the quantity of HSCs in human body. One might assume it 
to be, in a sense, proportionate to the average mass of a growing organism. Fig. H] demonstrates 
pooled data on the masses of human embryo, fetus, and body over time after impregnation. 
The data were approximated by a smooth curve 1, which is denoted M{t), then scaled linearly 
in log sea 1^ so as to move the first point to unit and last points to about 5 x 10^, known 
estimate of the number of HSCs in adults, curve 2. This dependence N(t) is used in further 
calculations. 

Time-dependence of the cell cycle time T(t) in humans may also be estimated from the 
age-dependence of the averaged human mass M{t). The duplication time at any moment of 
time follows from an obvious relation M{t) + M'i.{t)T^{t) = 2M{t). Then, 

As shown in Fig. [5], this function, computed from previously gained curve M{t), well ap- 
proximates points calculated from the empirical data on mass of the human organism. This 
calculation has been made based on a discrete analog of ([H]). 



^ Log-linear scaling is the simplest nonlinear scaling that transpose variable y io y' , values yi to y'^ and ?/2 
to j/2 so that this transposition is a linear one of a logarithm of y: y' ~ y[ exp [in (y/j/i) In (j/2/2/1) In"^ (2/2/2/1)] • 
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Figure 4: Human organism mass and the number of HSCs. Points: 
treatment of averaged data on masses of the human em- 
bryo, fetus, and body [137] . Lines: (1) approximation 
of the time dependence of the mass M(t), (2) scahng of 
the approximation curve M{t). 



Mass duphcation time Tj is not the cell cycle time T, though it directly relates to the 
latter. If cells of the human organism were dividing without death, T would be equal to Td, 
of course within a rough idealization in which a human organism is an ensemble of identical 
cells. However, one could assume that the regularity, which governs the exponential increase 
of the cell cycle duration, is approximately the same as for the mass duplication time. Then, a 
log-linear conversion would be a reasonable procedure to obtain T(t) from Td(t). The result of 
such conversion with 0.5 d and 10 d as initial and final values for T{t) is shown in Fig. [5l 

As an illustrative example only, we will take this function T(t) to study the degradation of 
HSCs, which are assumed to follow this dependence in early stage of the development of an 
organism. So, both functions N{t), number of cells in growing population, and T{t), cell cycle 
duration, are built on the tabulated dependence M{t) and fixed in the following calculations. 
This decreases the number of parameters, which is necessarily great in a model of such complex 
phenomena as cancer. On the one hand, even in this case we use initial and final values for 
and T. On the other hand, these values are natural for leukemic cells, and the semiquantitative 
results of the model are not sensitive to reasonable variations of these values. 

Two functions, and u{t), in ([7]) remain to be modelled to calculate probability of 
leukaemias. First, we note that even taking these functions constant gives rise to a peak in 
dP/dt, due to an interplay between the functions N{t) and T{t), curve 1 on Fig. El This 
might be associated with the observable growth of leukaemia incidence in early period of an 
organism's development. Position of the peak, some left of the empirical peak, depends on the 
tabulated curve N(t), precisely on the position on its rapid growth, and less sensitive to other 
parameters. This mismatch is not principal. Maximum peak value will be studied then. 

At constant u and u, the peak value is more than 10 times less than that of the empirical 
peak. This makes one think of another factors affecting the value of the local-in-time peak. 
The fact is that the peak value significantly depends on the character of u{t) and z/(t) time- 
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Figure 5: Time- dependence of the cell cycle time T. Treatment of 
averaged data on masses of the human embryo, fetus, 
and body [137] . 

dependences. 

According to medical knowledge about the immune system and DNA repair system their 
efficiencies grow in time in the early phase of an organism's development. On the whole, the 
immune system monotonously unfolds as cell differentiation occurs. The immune system gets 
stronger as a lot of previously unknown alien dangerous proteins become known. Apparently, 
in a similar way, the DNA repair system develops gradually as a lot of necessary protective 
proteins are expressed from specific genes. So the assumption, which is based on basic medical 
knowledge, is that the functions u{t) and z/(t) monotonously fall from an initial values to their 
final values, which means the efficiency of the immune and repair system grow correspondingly. 

The final values of u{t) and i>{t) have been already determined from the comparison with 
the empirical data. Fig. [21 within the frame of the model approach. Assuming that DNA 
repair system develops much faster than the immune system, we will take z/ as a constant and 
focus on the contribution to dP/dt from the immune system. If it were not so, synergistic 
effects might be great. It is convenient to model the function u{t) as a logistic curve u{t) = 
Ml + (mo - ui)/[l + {t/toY] with Ml = 1.6 X 10-^ uq = 7.6 x 10-^ s = 2, and to = 500 d, 
which gives a good fit. Fig. [6] curve 2. The changes in the peak value are shown with varying 
Uq, curves 1 and 3. Exponent s regulates the peak width. The peak value also changes with 
parameter to, but in a lesser extent. If to designate the maximum peak value p = max(dP/dt) 
and define a relative sensitivity to varying a parameter x by the formula S = {x/ p){dp/dx), 
then the sensitivity to parameters uq and a is close to unity, and that to to is about 0.37. 

Some discrepancy between the model and empirical data in the position of the peak could 
not be eliminated by varying model parameters in acceptable intervals. Most likely, the delay 
of empirical points vs calculated curve originates from the fact that it takes a few months in 
order that a single cancer cell, for which probability P is derived, would generate a detectable 
population of cancer cells. 

It would be interesting to examine how variations in the parameters, caused by chronic 
detrimental factors, like an enhanced free radical concentration, could affect the value of the 



10 



cu 
>, 



X3 
O 



10" 



10" 



10" 



10" 



O 4 




O 



0( 




10 



102 103 

t, days 



104 



Figure 6: Leukemia probability density dP/ dt calculated from (JTj): 
1, u = 2 X 10"^ and u = 1.6 x 10"^ as in Fig. El 2, 
uq = 7.6 X 10^''; 3, Uq = A X 10~^, explanation in the 
text. Age-specific leukaemia incidence rates in the USA, 
4, adapted from http://seer.cancer.gov 



peak in the probability density of cancer dP/dt. However, it is hardly possible. 

Excessive free radicals affect not only the current values of the parameters, but the course 
of their formation. It means it would be impossible to change the value of a parameter at a 
moment t without significant changing all the subsequent time dependence of that parameter. 
Only the rate of lesions could change directly with concentration of chronically acting genotoxic 
factors. The other parameters, z/, u, N, and T are systemic parameters. They describe the state 
of corresponding systems: the rate of DNA repair system, perfection of the immune system, 
the number of cells in population, and the averaged cell cycle duration. Each of these nonlinear 
functions depends on the prehistory of other systems as well. For example, the development of 
the immune system, the growth of the number of cells, and changes in the averaged cell cycle 
duration are evidently coupled into interacting processes. This makes it difficult to estimate 
quantitatively the consequences of a chronic extra content of free radicals, except a special case 
below. 

Clear qualitative results for the case of chronic factors may be obtained from analyzing 
equation (JTj). From ([71) follows that the time derivative of P, associated with age-specific 
cancer incidence rates, depend on several functions: 

dP 

—— = kp(a, z/, T)uNlT. 
dt 

It can be shown that the peak in dP/dt on Fig. [6l is formed mainly by two basic, opposite 
tendencies. Initial increase in dP/ dt is due to rapid exponential growth of the number of cells 
A^. This is opposed by decreasing u (growing perfection of the immune system) which finally 
prevails over the former tendency so that the peak is forming. The peak value strongly depends 
on whether the development of the immune system keeps pace with the growing number of cells. 
If it has not enough time or motivation to develop, a significant increase in cancer probability 
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may happen. Two strategies are evident to minimize the risk: one strategy is to promote 
formation of the immune system as soon as possible, the other one is to slow down the rapid 
growth of the mass of an organism. 

Unlike the effects of chronic factors, effects of short-term exposures to toxic factors may 
be easily calculated within the frames of this model. Consider again equation ([7]), from which 
follows 



Consider the value of P at a moment t' preceding the phase of the rapid growth of P, Fig. [HI 
As was shown, P{t') determines at which time the subsequent growth of P to dangerous values 
takes place. Time derivative dP/dt is directly proportional to any of the parameters capable 
of a sudden short-term change, a, z/, and u. Let x be one of them, and 5x be a short-term 
variation of it, localized at a moment t. Then the probability of a cancer to a moment t' and 
the increment 5P are equal to 



because time derivative of P practically remains unchanged during the time interval 5t within 
which 5x differs from zero. Evidently, 5P is the same for any t' > t + 6t. Quantitatively, 
increment 6P may be estimated if the data are known about the relative growth of a risk factor 
and about the age-specific rates of the cancer incidence. 

An important qualitative result: The long-term effect of a short-term exposure varies 
directly with the value of age-specific rate of cancer incidence on the moment of exposure. 
Leukaemia and many other cancers show enhanced incidence in early childhood. For example, 
based on the empirical data displayed on Fig. [6l one could conclude the following: medical 
x-ray procedures in the first few years are an order more hazardous, in terms of cancer prob- 
ability in old years, than those carried out in a later period. To our knowledge, it is the first 
mathematical observation of the correlation between elevated risk of these forms of cancer in 
old years and exposure to genotoxic factors precisely in the first years of life. 

Returning to chronic exposures, one could suggest a very rough estimate for the case when 
chronic variations 6x{t) of the parameters x{t) are small and change linearly with x, i.e. 6x oc 
x{t). In this case, the value 6x{t)/x does not depend on time. Then, conversions like in ffTUj) 
give rise to another relation for the increment 6P, which now depends on the time t' at which 
it is taken: 



i.e. 6P/P = 6x/x. Indeed, there are serious epidemiological and empirical data, as well as 
rationale, for a linear dose-response curve for radiation-induced cancers, particularly at low 
doses [95] . 

In the absence of any quantitative data on how small free radical variations could affect the 
different channels, that is, the rate of lesions and repairs, perfection of the immune system, 
and cell cycle durations, one might only assume, on the whole, in average, and as a very 
rough approximation, a linear dependence between corresponding relative values. That is, for 
example, a chronic 1%-change in the concentration of aggressive molecules causes 1%-change 
in each of the above channels. Contributions of small changes in different channels are summed 
up in P arithmetically, as was shown above. Moreover, since, on the whole, excessive free 
radicals do increase the rate of lesions, retard the DNA repair, impede the immune system to 
operate duly, and probably as a long-term systemic response to a danger stimulus, increase 




(9) 




(10) 
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the size of population and decrease cell cycle T, then, according to (E]) and (171), changes 
in all channels should be summed with the same plus sign. Consequently, there are grounds, 
though qualitative, to think that a small but chronically acting relative elevation of free radical 
concentration results in at least a few-fold greater relative growth of cancer probability P. 

Evidently, a 1%-increase in the probability of cancer in age 5-10 years will inevitably shift the 
curve P{t) up by the same 1%, even if no increase happens in the rate of genotoxic factors later 
and all the essential parameters remain the same. Fig. [3l It might seems to be an insignificant 
change. However, it means also that the curve P{t) moves left by about one month, which should 
be considered as a direct influence of P-level in early childhood on the lifetime. Therefore, 
analysis of possible physical reasons, which can increase probability of cancer in childhood, is 
well grounded. One of the possible reasons, not addressed in the literature until now, is magnetic 
nanoparticles chronically exposing nearby molecules to rather strong static MFs, which in turn 
promote formation of free radicals. 

3 Magnetic nanoparticles 

The agreement between the above model calculations and observable data is rather conditional, 
since many assumptions have been made on the character of key dependences. However, the 
model provides a possibility to study with greater certainty how different external factors might 
affect the growth of cancer probability. 

Among factors capable of significantly promoting cancers are well known ionizing and ul- 
traviolet radiations, some mutagenic chemicals related to tobacco use, environmental pollution, 
food contaminants, and some infectious agents. Other known risk factors are social factors: 
alcohol, obesity, unhealthy diet and physical inactivity, or factors intrinsic to an individual: 
age, sex, genotype and heredity. 

Last decades have shown that this list of risk factors might be widened. Recent epidemio- 
logical studies revealed a possible cancer risk of weak extremely-low-frequency electromagnetic 
fields generated mainly by power lines and various home electrical appliances. There is a corre- 
lation between the intensity of background magnetic fields and the rate of childhood leukaemia 
incidence. But the correlation is observed not always and another part of epidemiological stud- 
ies show no correlation. In spite of available controversy, the International Agency for Research 
on Cancer has classified magnetic fields as possible human carcinogens. A possible link between 
weak magnetic fields and cancer is hotly debated in several recent meta-analyses [2], [31], 11301 W2\ 
and reviews ^61], [871 127] . At the same time, the nature of the processes underlying this link re- 
mains unclear. In general, biological effects of electromagnetic fields, which do not appreciably 
heat biological tissues, look somewhat contradictory from the physical viewpoint, and there are 
no recognized physical mechanisms for such effects yet. 

A lot of hypothetical mechanisms have been suggested to explain biological effects of weak 
extremely-low-frequency MFs. A brief review of the mechanisms may be found in [9J and the 
detailed examination in [TU] . Most researchers often discuss the following hypothetical physical 
targets for MF action in magnetobiological phenomena: (i) iron-bearing magnetic nanoparticles 
growing in biological tissues, (ii) spin-correlated radical pairs, in some biochemical reactions, 
interacting with magnetic field by their spin magnetic moments, (iii) long-lived rotational states 
of some molecules inside protein structures, which interact with MF by their orbital magnetic 
moments. 

The basic problem is that the interaction energy of biologically active molecules and the 
MF at the geomagnetic level is very small. It is by many orders of magnitude smaller than the 
energy of thermal fluctuations fceT ^ 4 x 10~^^erg at physiological temperatures: <^ ^bT, 
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where n is the molecular magnetic moment and T is the absolute temperature, which is written 
in Roman to distinguish it from the cell cycle time T. Obviously, magnetic effects cannot exist 
here! It is the most heated argument raised by the opponents of the idea that weak MFs, 
on the order of the geomagnetic field and lesser, can affect organisms. At the same time, an 
ever-growing number of experimental evidences demonstrating biological effects of weak MFs 
require further attentive and careful theoretical study. 

Apparently, to explain observed biological effects of weak extremely-low-frequency magnetic 
fields one needs to equalize the inequality fiH <^ ^bT: either the magnetic moment fi of a 
suggested target in an organism should be sufficiently large, or the effective temperature T of 
the target should be sufficiently small. The former possibility is used in mechanisms based on 
magnetic nanoparticles |i53j found in tissues of many organisms including human brain tissues 
[HHllHl]. Magnetic energy of such particles may exceed k^T by many times and cause a biological 
response. In a recent work [8], as small MFs as of 200 nT were shown to significantly change 
the nanoparticles' dynamics. 

On the other hand, the above inequality relies on the implicit assumption that a target is 
in thermal equilibrium with the surrounding medium. Evidently, one could overcome the "kT 
problem" also having found possible targets whose effective temperature differs from that of the 
medium, so that the inequality would no longer be valid. Hence, no fundamental limitations 
would be placed on the possibility of observing biological effects of weak MFs that interact with 
such targets. 

Reactions involving free radical pairs give a clear example of the case where the inequality 
fails. Magnetic processes based on spin dynamics of the radicals develop so quickly that the 
thermodynamic equilibrium has no time to be established. This means spins move coherently 
and no temperature of spins exists within these small time intervals, for which the term "spin 
lifetime" is used. Another example where the inequality fails is the molecular gyroscope model 
[?]. In this model, a small biologically active molecule bound within a protein cavity is well 
isolated from the surrounding thermal perturbations. Interacting with the MF, it can coherently 
rotate for a long time. This makes the notion of temperature inapplicable to this type of 
molecular targets as well. 

As is seen, there are different ways to overcome the kT problem, at least at the conceptual 
level. Among them, the most certain are mechanisms involving magnetic nanoparticles as a 
primary target for MFs ^T\. A lot of studies demonstrate their presence in living tissues, on 
the one hand, and their involvement into biological reactions, on the other hand, [59] . However, 
possible role of magnetic nanoparticles in molecular processes of DNA degradation has not been 
discussed yet. 

Exposure of humans to airborne nanosized particles has increased dramatically over the 
last decades due to anthropogenic sources. Rapidly developing nanotechnology provides yet 
another source of potential contamination. Nanoparticles can penetrate across epithelial and 
endothelial cells into the blood and lymph circulation to reach bone marrow, lymph nodes, 
spleen, and heart. Entry to the central nervous system and ganglia via translocation along 
nerve cells is also possible [90j. Biological activity of nanoparticles stems from their specific 
surface properties and includes inflammatory, pro-, and antioxidant potential. Therefore, the 
uncontrolled presence of nanoparticles in the environment is now considered as a potential 
threat for human health within the frames of nanotoxicology |47j . 

Among others, magnetic nanoparticles possess very special properties because of their rel- 
atively large magnetic moment. In other words, such particles are small magnets that behave 
like a compass needle. On the one hand, magnetic nanoparticles can rotate in an external 
MF, thus exerting a pressure on biological tissues, and, on the other hand, they produce their 
own and relatively large MF, i.e. endogenous MF, which in turn may affect magnetosensitive 
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biochemical reactions. 

Magnetic nanoparticles may appear in an organism differently: (i) they can penetrate 
through the organism's surface as natural pollutants, especially as particles of iron oxides that 
abound in nature ; (ii) artificial magnetic nanoparticles may be introduced into an organism 
with certain aims and penetrate through cell membranes |1U2] ; (iii) magnetic nanoparticles can 
appear inside biological tissues in the course of natural process of biomineralization. 

Artificial magnetic nanoparticles are used in biology and medicine for different purposes. It 
is magnetic hyperthermia of tumors |117l HHl EOl [93], magnetic separation of proteins and cell 
sorting [106[ I104j . magnetically targeted drug delivery [16], delivery of plasmid DNA into the 
bacterial cells [IS], contrast enhancement in medical magnetic resonance imaging, and etc 

Below we discuss in what way magnetic nanoparticles might be involved into the DNA 
degradation process. Estimates will be made for the MPs produced by the particles and average 
increment of the free radical content due to the enhanced MPs. 

Magnetosomes have a magnetic moment and consequently they produce proper MP. This 
MP is not small, though it quickly drops with distance. An idealization is convenient, in which a 
magnetosome is a sphere of radius p with a point magnetic moment /j. in its center. In the point 
r = nr, the magnetic moment generates MP H = [3n(/xn) — /x] /r'^, where n is a unit vector in 
the direction of r. MP vector H has direction and magnitude. The proper MP direction near the 
magnetosome changes quite considerably with its oscillations and may change to the opposite 
direction with the oscillation frequency. However, inasmuch as we consider MP effects on the 
rates of the radical-pair reactions, the MP direction is practically of little significance. The 
averaging of molecular orientations of a radical pair (RP) usually occurs far faster than the MP 
direction alters. Therefore, the rates of these reactions are sensitive only to the MP magnitude. 
Note that an anisotropic contribution to the rate of RP-reactions is theoretically considered in 
[100] , but there are still no reliable experimental verifications for this hypothesis. Then, we will 
study the averaged values of a constant and alternating components of the absolute MP value, 
or MP magnitude. The absolute MP value follows from the above expression for a dipole field: 

H{r,e,ip) = ^Vl + 3cos2e , 

where fi = and r, 6, ip are the radius vector and azimuth and polar angles of the spherical 
reference system, and the vector fj, is directed along the Z axis. It is easy to find the MP of a 
fixed dipole averaged over the volume between two imaginary spheres of radii p and R: 

c'3p , R 
R^ — p-^ p 

where c' = 1 + ln(2 + y3)/(2\/3) ^ 1.38. Prom this we can readily find also the MP averaged 
over a sphere of radius r, if in the above formula for Hy we use the formal equalities p = r and 
R = r + dr, taking the limit as dr — >■ 0: H{r) = pc'/r^. Pig. [7] demonstrates the dependence 
of the magnetosome's average MP H, calculated by this formula, on the distance from the 
magnetosome surface. The fact by itself that the average proper MP of a magnetosome is 
many times greater than the geomagnetic field should result in some biological consequences. 

The MPs produced by magnetosomes vary in a wide range, but mainly they fall within the 
interval 1-200 mT. Do such constant fields cause biological effects in general? Many studies, 
in which biological systems were exposed to MPs of these magnitudes, demonstrated various 
biological effects. In particular, such MPs can increase the level of the DNA damages [521 11351 
I105j . deregulate the cell proliferative/apoptotic activity [28 | IIOH 136]. affect immunity [58] and 
gene expression [ 1271 13U] . Work reports on mutagenicity of MPs, tumorigenicity of MPs 
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Distance from magnetosome's surface, nm 

Figure 7: Proper MF of idealized magnetosomes of different radii. 

Its value is far exceeds the level of the geomagnetic field. 

was assumed in [H3], I125j . and MFs inhibitory activity with regard to apoptosis and tumor 
suppressor proteins was observed in |124l 1123] . 

It is not infrequently argued that these effects might explain the correlation between child- 
hood leukaemia incidence and the enhanced background extremely-low-frequency magnetic 
fields at places of residence |lU8t [UU] . 

A number of works report no effect of static MFs on organisms. It is difficult to recon- 
cile the above observations of different biological effects from static MFs with known reviews 
which claim absolute safety of medical magnetic resonance protocols |129[ 1110] . The magnetic 
resonance studies may really be safe due to their relatively short magnetic exposures. How- 
ever, this does not rule out the potential health hazard associated with chronic MF exposures 
and cumulative effects at the level of DNA mutations. Data are often contradictory: a part 
of studies report clear mutagenic, co-mutagenic, or toxic effects from tesla-range static MFs 
[Ha El EE], wWle others do not [ml[T86] . 

At least in part, this discrepancy may originate from backward influence of the genome 
on the magnetic effects: different strains of the same animals displayed significantly distinct 
sensitivity to co-carcinogenic or tumor-promoting effects of MF exposure [29] and different 
strains of the same bacteria showed clearly distinct frequency spectra of a magnetic effect [3]. 
Genes implicated in genetic instability syndromes were relevant in modulating the response of 
cells to an ac MF [76]. To add, the available epidemiological data are not sufficient to draw 
any conclusions about potential health effects of static MF exposure [30]. In this regard, a 
care should be taken: a negative result says either no effect takes place at all or the effect is 
less than the resolution of the experimental technique, a priori half by half. Usually, accuracy 
of measurements in biology is less than a few percent even for large statistical samples due to 
normally great inherent variability of biological objects. But a few percent effects are extremely 
important when considered with regard to human populations. 

Unlike constant fields, variable MFs induce electric currents in tissues, which may exceed 
natural biological currents, depending on the size and geometry of a biological system. It occurs 
when the frequency-amplitude product of an MF is greater than 10 HzxmT in order of value. 
For example, observation of a biological effect from 50 Hz MFs greater than 0.1 mT could not 
be followed by a statement that MFs of that magnitude entail the similar effects, because the 
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Figure 8: Pooled data on the relative biological effect from expo- 
sure to static magnetic field, adapted from [281 1123[ [52| 

1771 [miT^ . 

actual cause might be in the electric fields. 

Particularly heavy for interpretation are experiments, in which a biological system is exposed 
to an intermittent or pulsed MF. Within the short time intervals, when the MF quickly changes 
its magnitude, great electric current pulses appear in biological tissues, which undoubtedly 
causes significant electrophoretic-like biological effects. Here, the MF is not a prime cause of 
the effect, but just a method for inducing a great electric current in tissue, like in transcranial 
magnetic stimulation. 

For these reasons, all over the present article, we focus on the experiments, in which biolog- 
ical effects were produced by the constant MFs and extremely-low-frequency MFs lower than 
about 0.1 mT only. The following Table [T] summarizes the results of some recent experiments, 
sorted by the value of MFs used. 

As is seen, there are different effects of mT-range MFs on biological systems. Reviews 
[85l [22] discuss a variety of such effects with more details and in a wider range of MF inten- 
sities. Just a few experiments investigated the field dependence of magnetic effects. These 
experiments can provide information on the physical nature of the MF targets. Fig. [H] shows 
pooled MF dependences of the relative magnetic effect in different biological systems, in which 
the experimental data were normalized in a similar way. As a general motif of Fig. [HI one could 
assume that the relative biological effect is approximately proportional to the MF magnitude, 
reaching about 20-40% in 100 mT MFs. 

On the whole, the MF dependences of this type resemble those for the rate of magnetosen- 
sitive chemical reactions involving spin-correlated radical pairs. Radical pair mechanism would 
be a nice candidate to explain biological effects of weak geomagnetic-like MFs, unless some prin- 
cipal difficulties limited its sensitivity by maximum values of about 1%/mT for the majority 
of radical pairs |107l I118j . However, this limitation does not matter at all, regarding observa- 
tion of biological effects from MFs significantly greater than the geomagnetic field ~ 0.05 mT. 
There are experimental indications that radical pairs take part in primary interactions of mT- 
range MFs with biological systems [1031 11261 [TTl 1135] . Melatonin, known free radical scavenger, 
suppressed DNA damage induced by a static MF in [57j and by a relatively weak 60-Hz MF 
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Table 1: Effects of static MFs on biological processes 
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in [HE]- Other free radical scavengers blocked magnetic- field-induced DNA strand breaks in 
|69j : authors proposed that magnetic fields initiate an iron- dependent free radical generation 
process in cells, which can lead to genotoxic changes. The association between mT-range MFs 
and excessive production of free radicals is a robust hypothesis, since today no other equally 
elaborated idea exists of how such static MFs could affect biological processes. High priority is 
assigned by WHO to studies of a co-mutagenic effect concerning the carcinogenic potential of 
static magnetic fields with regard to RP reactions |132j . Thus, several reasons exist and many 
researchers are apt to think that biological effects from mT-range MFs, which vary roughly 
as the MF value, with a saturation in higher fields, originate from magnetochemical reactions 
with radical pairs, most likely highly reactive oxygen species [114^ 1115] . 

One of the possible sources of a chronic MF exposure is the enhanced content of magnetic 
nanoparticles in human tissues. Magnetic nanoparticles generate mT-range MFs and affect 
production of free radicals. The question is, how much? Let an idealized magnetosensitive 
reaction be AB AB ^ A + B, where the intermediate AB is a spin-correlated RP in a "cage" 
formed by the molecules of the surrounding viscous medium. 

One of the results of magnetochemistry is that the relative rate of recombination AB AB 
may change between 1 and 1/4, in some units, depending on the state of spin evolutions in the 
RP. An MF may change the rate just within this limits. Usually, magnetic field dependences in 
magnetochemistry are complex dependences that may contain minima and maxima. However, 
as a whole, the magnetic effects grow with MF. As a very rough approximation, we assume 
that the rate of free radical formation K, which is bound to the recombination rate, depends 
on the value of magnetic field H in the following way 

K = Ko{l + 3 „^ „ ] , (111 



Ho + H 

where Kq is a minimum rate at if = and Hq is a characteristic MF, which lies between 10 mT 
and 1 T for most RP reactions. The approximation is justified also due to a lot of different RP 
reactions in organisms, which are poorly identifiable in experiments on MF effects. We note, 
that the above formula is an optimistic one: usually, just 1-2%-relative changes are observed 
even in strong MFs. However, looking at Fig. [HI one could assume that a few percent changes 
in free radical content should lead to about ten fold higher alterations in biological observables, 
probably due to an interplay between nonlinear space- and time-scaled biochemical reactions. 

Averaged value of K over a volume of radius R surrounding a magnetosome of radius p 
is of interest. In calculations, the averaged MF of a magnetosome H ^ ix\^/r^ is taken for 
H and the magnetic moment /i = vJ, where v and J are the magnetosome's volume and the 
saturation magnetization of magnetite Fe304. Averaged value of K, in Kq units, is then as 
follows 

^^~^ + a(iiVp^-l) '"^T7^- -3ff„/(4.v/2J). (12) 

This quantity is plotted on Fig. [9] as a function of R/p at different values of a. As is seen, the 
average rate of the reaction markedly exceeds the level 1, which corresponds to the absence of 
the MF, only within the space limited by about ten-fold magnetosome radius. 

It is known that for the majority of RP reactions, at the RP lifetime r of the order of 10~^ s, 
the characteristic field Hq is about 100 mT in the order of value. Correspondingly, the value 
of the parameter a in ( fT2i) equals about 0.1. For estimations, an inverse proportion between 
r and Hq may be assumed, so that relatively long lifetime r ~ 1 /is means that a equals 
10~^. However, so large lifetime could be an exceptional case: about 20 ps for the flavin- 
tryptophan [211 EH] and 2-10 ps for flavin-tyrosyl |131j . Therefore, it would be reasonable to 
assume a ~ 0.1, and thus the rate of most RP reactions might be shifted significantly only near 
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Figure 9: The average values of the rate of an ideahzed RP reac- 
tion in a volume of radius R surrounding a magnetosome 
of radius p. 



magnetosomes, within the distance R ^ p from magnetosomes' surface. At the same time, the 
presence of magnetic nanoparticles means a chronic magnetic exposure. In this case, as was 
shown above, even a few percent changes in the rate of RPs formation might be biologically 
significant due to accumulation of their contribution to DNA degradation. 5-% changes in K 
take place over as large volume as of the radius R ~ 10-20 p for the straightforward case 
a = 0.1 and R ~ 100-200 p for a = 10~^. This unequivocally indicates that if a cell contains a 
magnetic particle, even as small as of a few nanometers, a significant part of the cell is covered 
by the MF of that particle, so that the rate of RP formation grows by about a few percent. 

Probably, there is a hidden association between magnetic nanoparticles, which are composed 
of iron oxides mainly, and iron ions in cytoplasm that were shown to significantly increase the 
number of rat lymphocytes with DNA damage when administered simultaneously with the 
exposure to static MF of 7 mT [56]. The role of iron ions in free radical production is discussed 
in [69]. 

There are some proteins, like ferritin accumulating several thousands of iron atoms, which 
are necessary for biogenic magnetite production from ferrihydrite they contain [72| HOj. Iron 
compounds related to magnetite and ferritin were observed in a tumor tissue in [H]. Nanoscopic 
magnetic crystals in ferritin cores may be a byproduct of biochemical processes utilizing Fe 
ions [98]. One might hypothesize that another yet unidentified proteins control and impede 
magnetite production, since its presence in extra amounts is harmful or toxic. Then, an equi- 
librium size of the magnetite byproducts should exist, since counteracting stimulus are to be 
proportionate to the surface of magnetite substance, which grows with size slower than mass. 
Probably, the size is very small and concentration of magnetic particles is extremely low under 
normal conditions, so that no biological consequences are possible. However, the dynamical 
equilibrium may be disturbed by the absence or inactivity of the specific proteins. As well, 
the equilibrium may vary in different organisms by genetic reasons, resulting in a chronically 
increased or decreased mean free radical concentration and thus in a higher or lower rate of 
genetic degradation and different effects on cancer probability. 

Time-course of the formation of extracellular biogenic magnetite was analyzed by quantum 
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magnetic measurements in the culture of Geobacter metallireducens |128j . Authors of [109] . 
using cryo-electron tomography that allows 3D reconstruction at rather high resolution, moni- 
tored the process of magnetosome formation in M. gryphiswaldense. These studies have shown 
that transformation of ferrihydrite to magnetite is a complex many-stage process that can lead 
to formation of 5-10-nm crystallites within 30 min. Larger crystals grow during several weeks. 

Magnetic nanoparticles, depending on their size and substance, are in a multidomain, single- 
domain (SD), or superparamagnetic (SP) state. Biological formation of multidomain particles 
is seldom occurred. SD magnetite particles are of about 15-100 nm in size. They have received 
most attention in scientific literature as it is the most obvious and relatively easy to explain 
how they move mechanically, how external MFs act on their motion, and how this gives rise to 
activation of mechanoreceptors. The magnetic moment of an SD particle is rigidly bound to 
it, so that the MF interacting with the magnetic moment exerts a torque on the particle and 
mechanically rotates it. 

Unlike SD particles, SP particles have their magnetic moments mostly unbound from their 
geometry. The magnetic moment of an SP particle can be switching between a few "easy" 
directions, which are determined by minima of the total magnetic energy. The latter includes 
the magnetic anisotropy energy of the bulk substance and demagnetizing energy that depends on 
the particle shape. Thermal perturbations make the particle to chaotically change its magnetic 
moment orientation, and no magnetic moment appears on average over time. An external MF 
brings the additional magnetic energy to the particle and makes it to preferably orient in only 
one of the easy directions. The time-averaged magnetic moment appears, however in weak 
MFs, it is very small as compared to the instant magnetic moment, and its magnetic energy is 
many orders of value lower than /cbT. For this reason, no possible effects of SP particles have 
been studied thus far, except a few models where macroscopic conglomerates of such particles 
were considered from the viewpoint of their magnetic interaction with each other to explain 
magnetic navigation in some birds [20j . 

In this article, an important role of single SP magnetic nanoparticles is emphasized. There 
are different ways in which they could magnetically influence biologically significant processes, 
in particular, the direct influence of the proper MF of magnetosomes on free radical formation, 
considered above. Another way is connected with the ability of SP nanoparticles to switch their 
magnetic moments time to time, being driven by random thermal forces. Frequency spectrum 
of this change-over depends mainly on the energy of the potential barriers, or maxima of the 
total magnetic energy. Since the switching is conjugated with a quick magnetic flux change, an 
electric eddy current, localized near the particle, is induced by this switching. On the one hand, 
this electric current should affect nearby chemical processes. On the other hand, this presents 
a dissipative factor for the switching process. Therefore, the magnetic switching, a non-inertial 
process in principle, is described by a dissipative dynamics with a nonlinear potential function. 
This is very similar to the stochastic dynamics of SD magnetosomes considered in [8] . Different 
dynamic conditions are possible here, including stochastic resonance in external ac MFs and 
switching frequency control by static MFs. 

The density of magnetosomes in the human brain is more than 5 x 10^, and in meninges 
more than 10® crystals per gram [03]. In fact, about 90% of the particles measured in this 
work were 10-70nm in size, and 10% were 90-200 nm. Subsequent studies have shown that 
the concentration of magnetite/maghemite in human tissues varies from tens to hundreds ng/g 
[H] and equals about 50 ng/g on average in the human brain |112j . Magnetite levels may be 
even higher when considering superparamagnetic particles, diseased tissue and age [211 HE] • No 
data are known about possible content of SP nanoparticles in human blood or marrow. Their 
possible contribution to magnetically measurable quantities could be masked by relatively large 
amounts of ferritin [15] and diamagnetic and paramagnetic fraction of blood, which requires 
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development of special methods for single SP nanoparticles detection. 

The SP particles present a new and highly interesting object from the viewpoint both of 
quantum physics and of biology, particularly biology of cancer. Specific molecular mechanisms 
based on SP nanoparticle dynamics will be studied elsewhere. 

4 Discussion 

Enhanced cancer probability in the early period of life results from rapid growth in the number 
of cells in a population while the repair and immune systems are not yet ready to repair DNA 
damage effectively and remove the cells with mutant DNA molecules. 

To be protected from early cancer, it is very important for an organism to form the DNA 
repair system and immune system as soon as possible, since even a small delay may result in 
a significant growth of risk. Therefore, exogenous factors may be dangerous even when they, 
per se, do not cause DNA lesions or mutations. If exogenous factors retard the formation of 
protective systems, they promote enhanced cancer risk. 

Formation of the immune system may go faster under training, i.e., during some slight viral 
diseases or other infections. Paradoxically, developing organisms need to be sick from time 
to time in order to be healthy on the whole. There is recent epidemiological evidence that 
any activity outside the home during the first year of life, especially the first three months, is 
associated with a decreased incidence of acute lymphoblastic leukaemia, by 40-50% [BB]; the 
most likely explanation being a protective effect from exposure to common infections. It is 
interesting that the mathematical model demonstrates this effect. 

The other non-evident conclusion is that the risk of early leukaemias might be reduced by 
slowing down the division of cells of a growing organism in the first year(s) after impregnation in 
order to give more time to the immune system to develop. It is probable that a rough life for a 
woman during pregnancy and for her child after birth could promote an advanced development 
of immunity and decrease the probability of cancer in the first years. Children conceived in 
lean years may show a reduced cancer incidence rate in early age. 

The DNA degradation model described here is a probabilistic model. Though valid with 
respect to cancer risk assessment for individuals, it could be verified only by epidemiological 
studies. However, the situation may change in the future. Assuming that the probability of 
cancer P may be determined not only for a human population, as P{t) in Fig. [3], but for 
a separate human organism as well, one might speculate that there is an oncological age. 
Formally, it is the age t where P{t) equals to the probability of cancer for that human being. 
Oncological age, unlike biological age and, moreover, physical age, can be significantly affected 
by long-term detrimental environmental conditions and short-term risk factors. Oncological 
age, very similar to Balzac's magic piece of shagreen, shows the probability of cancer for this 
human in terms of age. Most likely, the future success of nanotechnology will make it possible 
to straighten DNA molecules and read out genetic information quickly by physics methods. 
Then, oncological age will become a truly measurable quantity. 

The model might be further corrected in many aspects. An obvious one is that the model 
has no specificity with regard to different genes. One more aspect is that there are different 
classes of DNA damage /mutations and they have different rates. Another aspect is that the 
model studies a population of idealized identical cells, while cells with parameters, distributed 
over some ranges, would be more realistic. For example, it would be right to assume that 
only some of the cells contain magnetic nanoparticles and then have a significantly increased 
rate of DNA degradation, rather than that all cells have a moderately increased rate. Yet 
another aspect is that there are "hot and cold spots" for spontaneous mutations, i.e. DNA 
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sites that have a higher or lower rate of mutations than expected from a normal distribution. 
Apparently, the list might be continued. However, even in its current simple form, the model 
is effective. Basic clinical observations, like the exponentially growing probability of cancers in 
middle age and an enhancement of the probability for some cancers in early childhood, have a 
clear explanation within a modeling context. 

Some of the results of the model are probably new and unobvious as mathematical corollar- 
ies: there is a synergistic effect of harmful factors, like free radicals, acting through the different 
biological/biochemical pathways or channels of DNA degradation, and there are unequal long- 
term consequences from the same short-term detrimental exposure applied in different phases 
of an organism's development. 

A possible source of an increased level of free radicals is exposure to relatively strong MFs 
from magnetic nanoparticles in human tissues. If a cell contains a magnetic particle, even 
as small as a few nanometers across, a significant part of the cell is covered by the MF of 
that particle, so that the rate of free radical formation may grow by a few percent. Magnetic 
nanoparticles, chronically increasing local concentration of free radicals around themselves in 
cells, may affect cancer related processes through all the channels: increasing rate of lesions 
to DNA molecules, decreasing rate of DNA repair, impeding immune responses and apoptosis, 
causing acceleration of an organism's development and shortening cell cycles. Since, in general, 
biological effects of MFs are known not to be strong, it is reasonable to assume that a most 
sensitive system is a more probable target for MFs in living matter. Within the model context, 
it is the immune system. 

External ac MFs may change the dynamics of magnetic nanoparticles in different ways. 
By facilitating the switching of the particle magnetic moment and so the rate of electric pulses 
induced, ac MFs can increase the rate of free radical formation. By shaking magnetic moments, 
ac MFs can impede adaptation of cell structures to their presence. It can also shift the spectrum 
of the endogenous ac MFs, and unbalance time-ordered cycles of RP biochemical reactions. 
These specific mechanisms are hypothetical and require further estimates. 

How to test whether SP nanoparticles contribute to genome degradation? The most direct 
way would be to observe a difference in the amount of the particles in tissues of cancer patients 
and healthy people. Another approach to identifying this specific risk factor is to search for an 
association between the amount of the particles and the age at diagnosis of a leukaemia. Some 
indirect indications may come from epidemiological studies. On the one hand, stochastic reso- 
nance of magnetosomes in ac MFs [8] shows non-trivial static MF dependence of the induced 
biological effects. On the other hand, an increasing trend of childhood leukaemia were associ- 
ated with certain combinations of household static and power-frequency MFs [13]. Laboratory 
studies could also be useful: SP nanoparticles may be artificially introduced into cell cultures 
or whole organisms, which could result in cancer development in chronic experiments. 

5 Conclusion 

In summary, a mathematical model has been developed that explains basic oncological reg- 
ularities and offers unobvious strategies for minimizing risks of some cancers. It is proposed 
and substantiated that an enhanced level of acute lymphoblastic leukaemia in early childhood 
may originate from magnetic nanoparticles located in hematopoietic stem cells. The magnetic 
nanoparticles may have a natural biogenic origin from ferritin or they may appear in the cells 
due to exogenous contamination by widely spread iron oxide nanoparticles. Alternating mag- 
netic fields interacting with these particles change their contribution to free radical biochemical 
processes, responsible for the development of the immune system. Excessive free radicals impede 
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development of this system and delay its maturation, thus resulting in an enhanced leukaemia 
incidence. 
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Addendum. Numerical values of the model parameters 

Human genes vary enormously in size and exon content. Exons are regions of genes that code for 
main portion of proteins. There is an inverse relationship between gene length and percentage 
of exon content, so that a great variety of different genes have exon lengths between 0.2 and 
20 kb [120J. Apparently, only mutations to exons have detectable phenotypic effects, while 
others do not alter protein functions. Due to these facts, the parameter k, the number of bits 
in a gene is taken to be /c = 2 x 10^, i.e. in the middle of the range 0.2-20 kb. 

Mutations are induced in DNA by exposure to a variety of mutagenic factors in both in- 
tracellular and external environment. Spontaneous errors in DNA replication and repair are 
the main source of mutations. For higher eukaryotes, the rate of spontaneous mutation per 
generation is about one mutation per gamete |120j . The size of the human DNA sequence is 
near 3 x 10^ bp. So, the rate of mutations per bp, or bit, is about p ~ 3 x 10~^°. Since the 
coding DNA of an average human gene is 0.2-20 kb, mutations occur spontaneously with an 
average rate from about 6 x 10~^ to 6 x 10~^ per gene per cell division, which is in accordance 
with known estimates of the rate of spontaneous mutation [251 EH [26]. In the present model, 
the probability of mutation per bit per division p ~ au that follows from ([Tj). Since z/ is a fitting 
parameter, found to be about 10~^ d to fit data on Fig. [21 parameter a, the rate of lesions per 
bit per cell division, has been taken 10~^. 

HSCs enter the cell cycle with different rates, since in vivo they are usually in different cell 
cycle phases. In mice, approximately 8% of HSCs asynchronously entered the cell cycle per day 
and about 50% in 6 days [12] . Cell cycle parameters of murine bone marrow cells were estimated 
in m based on the comparison between a cell proliferation model [73] and experimental data 
from cell tracking experiments. Cell cycle duration was found to be about 0.3 d. In our work, 
the cycle duration of HSCs is a function of time within the limits 0.5 to 10 d, obtained by a 
conversion from the tabulated time dependence of the mass of the human organism. 

Reparation time r describes the efficiency of the DNA repair system. The repair efficiency 
is assumed to depend on accumulated mutations. At the initial relatively long period of a 
slow growth of P, time r is practically equal to its initial value u. Apparently, there are a 
wide spectrum of times needed to repair DNA damages. Normal cells react to an unrepaired 
DNA damage by retarding the cell cycle at a checkpoint until the damage is repaired, and 
triggering apoptosis when the damage is lethal. Retardation may continue from days to years. 
Proliferation kinetics of HSCs isolated from mice has been studied in [TF]: about 8% of cells 
asynchronously entered the cell cycle per day, 75% of cells were quiescent at any one time, and 
all cells were recruited into cycle regularly such that 99% of them divided on average every 57 
days. The values between 20 and 50 days were used in [75] in modelling HSC kinetics. 

It is unclear to which extent one could identify retardation time and repair time. Complex 
"digital" interaction between repair process, cell cycle arrest and apoptosis has been described 
in [81J. Two character times, fast and slow, of the double-strand break repair process were 
suggested in |119] : most of the breaks are rejoined quickly for 15 min and the rest are rejoined 
very slowly, for 10 to 15 h repair half-time. In the present model, cells enter the next cycle 
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without delay, and repair time is a separate intracellular process. The value of repair time u 
was taken on the order of 0.001 d, i.e. about a minute, to get probability curves fitted well 
to epidemiological data. Fig. O This means, in this model, parameters v and r should not 
be deemed precisely as repair time, it is rather parameters describing overall efficiency of the 
repair system. 
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